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Abstract 

We demonstrate that the analytic calculation of the 1RSB break point parameter in a paper 
by de Oliveira and Fontanari|l) is erroneous, due to the omission of a higher order term in 
a lengthy perturbative calculation, and provide a refinement of the accompanying numerical 
results. 

In 1999, de Oliveira and Fontanari (OF) studied the one-step replica symmetry breaking 
(1RSB) of a glass of Ising spins with a quenched random p-spin interaction of infinite range 
in a field Q. The Hamiltonian is given by 

U= Jii-ipVh ■••(?%„ -h^2<7i (1) 

«i<«2-"<*p i 

where the are independent Gaussian variables with zero mean and variance p\J 2 /2N P ~ 1 . 

They found that for fields h less than a critical value h c the transition was discontinuous (D1RSB), 
while for h > h c it is continuous (C1RSB). In §3.1 of that paper, they give certain results on the 
C1RSB line, including an expression for the 1RSB break point quantity x (m in their notation). 
We demonstrate an error in this calculation of x. Moreover, they present numerical results within 
the 1RSB phase which are (for reasons we shall explain below) inaccurate near the C1RSB line. 
We present a discussion and refined results. 

The self-consistent equation for the RS phase is|| 

q = T(2) (2) 

where 

T(n) = J -jL e- z2/2 tanh" \pf3 2 J 2 q^ z + j3h\ . (3) 

This solution becomes unstable against small replica symmetry breaking fluctuations on the 
Almeida-Thouless line, given by[Q 



where 



kS(4) = 1 (4) 



k=\p{p-l)0 2 J 2 <f-\ (5) 
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where 



G = sJlp^J^- 1 z + ^JlppJ^- 1 - eft- 1 ) z 1+ !3h. 



= (7c) 



(8) 



The C1RSB line is defined by <Zo = <Zi- In this case, G does not depend on z\, and the integrals 
over this variable are trivial. We observe that (|7a|) and (Q) reduce to the same equation, namely 
(||): on this line, the solution coincides with the RS, as we would expect. We obtain a second piece 
of information by subtracting these equations and performing a series expansion in the quantity 
e = qi — qo: (|7a|) and ( fTq ) both become 



and the difference (|7b| ) — (7a) becomes 



q = T(2) + 0(e) 



kS{A)e + 0{e) 2 



(9) 



(10) 



So (Q) is also satisfied on the C1RSB line, that is the transition coincides with the onset of 
instability in the RS solution, again as we would expect. 

On the C1RSB line, J7c] ) is trivially solved. We obtain further information by a series expansion. 
To first order we get 



l -kqe= l -kT{2)e + 0{ef. 



(11) 



This simply tells us that (j^) holds, which we already knew. We therefore eliminate this first order 
term by subtracting ( [7a] ) multiplied by ke/2 from (^) to obtain a new equation. To second order 
we get, rearranging slightly, 



-[g-T(2)] + [l-*S(4)] \e* + 0{e) 



0. 



(12) 



This tells us that (0) holds, which again we already knew. We therefore eliminate these second 
order terms by subtracting ( [7a] ) multiplied by (p — 2)ke 2 /Aq and [(^) — (|7a|)] multiplied by ke/A 
to obtain another new equation. To third order we get 



2Aq 



-C + AqkS(A) - 6qkS(6) + 2qkS(6)(l ~ x) 



e 3 + 0{ef = 



where 



G = 2(p - 1) - 2{p - 3)-^ - 3fcS(4) 



(13) 



(14) 
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We solve this to obtain an expression for x near the C1RSB line: 



6gfcS(6)-4gfcS(4)-(p-2)Cfc- 1 
1 ~ 2; - ZkqW) +0(e) (15) 

_ 6gfcS(6)-[4gfc+(p-2)]S(4) | 
2kqS(6) 

where we have used (^) and (|l^) to simplify our expression. This differs from (36) and surrounding 
equations of OF. We note that if one erroneously neglects the O(e) terms of ([7a]) when subtracting 
that equation multiplied by (p — 2)fce 2 /4g in the above process (using only the leading order 
equation (|^) instead) one obtains an incorrect form of the 0(e) 3 equation ( |l3| ) which gives exactly 
the form of OF. 

Since we have shown the analytic expression of OF for x near the C1RSB line to be incorrect, 
we must question the accuracy of their numerical solutions of the 1RSB equations (^) in that 
region, as the latter appeared to corroborate the former. The determination of x is indeed rather 
delicate, as the x-dependence of these equations is very weak, for reasons that are clear from the 
above analysis: it appears in a term 0(qi — go) 2 smaller than the leading order, and close to 
C1RSB, (gi - g ) < 1 by definition. 

We adopt an approach designed to avoid this problem. Rather than solving the equations as 
given, we choose a judicious linear combination which does not possess the same flatness. We 
know from above that subtracting ( |7a| ) multiplied by fc(gi — go)/2 from ( |7c| ) eliminates the leading 
order of the latter, leaving an equation where the x-dependence is suppressed only by a factor 



0(qi — go); and that further subtracting (7a) multiplied by (p— 2)fc(gi — go) 2 /4g and [(7b) — (|7a|)] 
multiplied by k(q\ — go)/4 eliminates the next order, leaving an equation whose leading order is 
linear in x. We find it most efficient to use the second of these very close to C1RSB (where the 
problem is worst) and the first elsewhere. We use a modified form of Newton-Raphson to find the 
roots, and do the numerical integration using Gauss-Hermitc quadrature. 

The figure shows the solution at p — 3 and h/J = 1 as a function of temperature. This is 
equivalent to the solid line in figure 3 of OF. The main line shows the numerical solution for x. The 
two diamonds show the predictions for x on the C1RSB line, the lower using ( |l6| ) and the upper 
using the equivalent expression of OF. The inset shows an enlargement of the region around the 
transition, with rectangles for the numerical predictions for x and a diamond for our perturbative 
result. 
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